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Abstract 



\ The perturbative QCD fragmentation functions for a heavy quark to fragment into 



heavy- light mesons are studied in the heavy-quark limit. The fragmentation functions 
for S-wave pseudoscalar and vector mesons are calculated to next-to-leading order in 
the heavy-quark mass expansion using the methods of heavy-quark effective theory. 
The results agree with the — > oo limit of the perturbative QCD fragmentation 
functions for b into B c and B*. We discuss the application of the perturbative QCD 
fragmentation functions as a model for the fragmentation of heavy quarks into heavy- 
light mesons. Using this model, we predict the fraction Py of heavy-light mesons that 
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are produced in the vector meson state as functions of the longitudinal momentum 
fraction z and the transverse momentum relative to the jet axis. The fraction Py is 
predicted to vary from about 1/2 at small z to almost 3/4 near z = 1. 
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I. Introduction 



Heavy-quark spin-flavor symmetries are very useful for understanding the properties of 
hadrons containing a single heavy quark in kinematic regimes where nonperturbative aspects 
of the strong interaction are dominant. These symmetries arise from the fact that the charm, 
bottom, and top quarks are much heavier than Aqcd- The symmetry is exact in the limit 
of infinite quark mass, and corrections can be systematically organized into an expansion 
in powers of Aqcd/^q using heavy-quark effective theory (HQET). There has been much 
progress on the applications of heavy-quark symmetries and HQET to the spectroscopy, and 
to both exclusive and inclusive decays, of charm and bottom hadrons |T|. 

It has recently been pointed out by Jaffe and Randall []2| that HQET can also be applied 
to the fragmentation of a heavy quark into hadrons containing a single heavy quark. They 
showed that when the fragmentation function is expressed in terms of an appropriate scaling 
variable, it has a well-defined heavy-quark mass expansion. Specifically, they showed that the 
fragmentation function Dq^h(z) at the heavy-quark mass scale has a systematic expansion 
in inverse powers of tuq when expressed as a function of the scaling variable 

rz 

where r = (m# — mQ)/m,H, is the mass of the heavy hadron, and z is its longitudinal 
momentum fraction relative to the fragmenting heavy quark. In the case of a heavy-light 
meson, r can be interpreted as the ratio of the constituent mass of the light quark to the 
meson mass. For the pseudoscalar meson P and vector meson V of the same S-wave multiplet 
( 1 S'o, 3 5'i), the fragmentation functions at the scale tuq have heavy-quark mass expansions 
of the form 

D Q _> P {z) = a ^- + b(y)+0(r), (2) 
r 

D Q ^v{z) = a ^- + b*(y) + 0(r), (3) 
where a*(y) = 3a(y). By heavy-quark spin symmetry, the leading terms differ by a spin 



factor of 3, while spin splittings first appear at next-to-leading order in the functions b(y) 
and b*(y). 

It was also realized recently that the fragmentation functions for mesons containing a 
heavy quark and a heavy antiquark can be computed using perturbative quantum chromo- 
dynamics (PQCD) [[3], [|, |!|. The fragmentation functions for a b to split into the S-wave be 
mesons B c and B* were calculated to leading order in a s in Ref. ||. These fragmentation 
functions have been used to predict the production rates of the B c meson at LEP and at 
the Tevatron |7], §]. The general analysis of Jaffe and Randall must certainly apply to per- 
turbative QCD fragmentation functions in the limit where the mass of the heavier quark is 
taken to infinity. It was verified explicitly in Ref. |]] that the PQCD fragmentation functions 
Di_> Bc (z) and Di_, B *(z) reduce to the forms (0) and ([|) with r = m c /(mb + m c ) in the limit 



Since the PQCD fragmentation functions are consistent with heavy-quark symmetry, 
they can be used as models for the fragmentation of heavy quarks into heavy-light mesons. 
In this paper, we show how the leading and next-to- leading terms in the I/tuq expansions 
can be calculated directly from HQET. We then discuss the use of the PQCD fragmentation 
functions as a phenomenological model for the fragmentation of charm and bottom quarks 
into heavy-light mesons. As an application of this model, we consider the fraction Py of 
heavy-light mesons that are produced in the vector-meson state. 

II. PQCD Fragmentation Functions from HQET 

The HQET Lagrangian, including the leading and the 1 / mq terms, is given by 



oo. 



£ = h v liv ■ D + 




2m Q 



1 




V 



(4) 
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where 

Ci = 1, 

ft - 3(^r 8/<33 ~ 2 ""-2. < s > 

g,(ri ^-9/(33-2^) 



These coefficients are all equal to 1 at the heavy-quark mass scale /i = tuq. The term propor- 
tional to C*2 can be omitted in calculating physical quantities, because it can be eliminated 
using a field redefinition involving the equation of motion (v ■ D)h v = from the leading 
term in the Lagrangian. Our method for calculating the fragmentation function involves a 
heavy quark which is off-shell by an amount at least of order iriQm q . To demonstrate that 
the C*2 term can still be omitted in this case, we keep it in our calculation throughout, and 
show that it cancels between the vertex and propagator corrections. 

In our calculation, we need the Feynman rules derived from the HQET Lagrangian for (i) 
the heavy-quark propagator, including l/rriQ corrections, (ii) the heavy-quark-gluon vertex, 
including 1/mg corrections, and (iii) the propagator for the small component of the Dirac 
field of the heavy quark. This last Feynman rule is needed in our calculation because 
the fragmenting heavy-quark is off its mass shell. The Feynman rule for a heavy-quark 
propagator is 

vk + ^k 2 -^(vk) 2 2 ' 1 ' 

2m Q 2mQ V / 

where k is the residual 4-momentum of the heavy quark. The QQg vertex is 

- ig s T a (V + -— !-(Aa + k 2 y - ^-v ■ (h + fc 2 K + i^-^Qu) , (7) 
V 2m Q 2m Q 2m Q J 

where k\ and k 2 are the residual 4-momenta of the incoming and outgoing quarks and 
q = k 2 — k\ is the momentum of the gluon. The Feynman rule for the propagator of the 
small component of the Dirac field of the heavy-quark is 

V^k — q ») — ■ (8) 

To calculate the fragmentation functions, we follow the method introduced in Refs. 
and H and applied in Ref. [|] to the fragmentation processes b — >• B c and b —>■ B*. We 



denote the pseudoscalar and vector Qq mesons by P and V, respectively. Here Q is the 
heavy quark, and q is the light antiquark. We calculate the cross section for producing a 
Qq meson plus a light quark q with total 4-momentum K^, divide it by the cross section for 
producing an on-shell Q with the same 3-momentum K, and take the limit K Q — > oo. The 
fragmentation function is 

D(z) i f dse ( s .Ml.J±) lim (9) 



167T 2 J \ Z 1-zJ K -+oo \Mo\ 

where M = tjiq +m q is the mass of the meson in the nonrelativistic approximation, s = K 2 , 
Ai is the matrix element for producing P + q or V + q, and M.q is the matrix element for 
producing an on-shell Q. The calculation can be greatly simplified by using the axial gauge 



with the gauge parameter n M = (1, 0, 0, —1) in the frame where = (Kq, 0, 0, y K 2 — s). In 
this gauge, we need only consider the production of the Qq meson plus q through a virtual Q 
of momentum K 11 . The part of the matrix element M. that involves production of the virtual 
Q can be treated as an unknown Dirac spinor V . In the limit K — ► oo, the same spinor factor 
T appears in the matrix element M.q = Tu(K) for an on-shell Q. The Feynman diagram for 
Q* — > Qq + q is shown in Fig. |I| The usual projection of the Qq onto a nonrelativistic 1 5o 
bound state reduces in the heavy-quark limit to the Feynman rule 

Qq ^rf^fe-f2> (10) 



where R(0) is the radial wavefunction at the origin for the meson and is its 4-velocity. 
For the 3 S\ state, the projection is the same except that 7 5 is replaced by j-, where e M is 
the polarization 4-vector for the vector meson V. The rest of the amplitude corresponding 
to Fig. |l] is obtained by using the ordinary QCD Feynman rules for the light quark spinor 
and the qqg vertex and HQET Feynman rules for the heavy-quark propagator and the QQg 
vertex. 

The amplitude M. for producing the 1 Sq state, including 1/mq corrections in the heavy- 
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quark propagator and vertex, is 



8 v /7ra s i?(0) M 2 y/M 1 1 / n^K + k^ 



i-M — _ (n _ 2Wl , Ciniq Ca , \9nv 

mg 2mQ 



m q (s - m 2 Q ) 2 l + _ JZl v . fc ^ n ■ k 



x n(p')7V(l + f) U + -^k" -^(v- k)v v + -^(7^ - W)) ^ T (ll) 
\ 2m Q 2m Q Am Q J 2 

where k = m q v + p' is the momentum of the virtual gluon and also the residual momentum 
of the fragmenting heavy quark: K = rrigv + k. Note that the term proportional to n u in 
the numerator of the axial-gauge propagator for the gluon vanishes after contracting with 
the Dirac factor. For the vector meson state, the 7 5 in the above equation is replaced by 

We are interested only in the sum of the first two terms a(y)/r + b(y) in the heavy-quark 
mass expansion, where r = m ? /(mg + m q ). We calculate separately the contributions to 
the fragmentation functions from the leading terms in the HQET Feynman rules, from the 
l / mq corrections from the propagator, and from the l / itiq corrections from the vertex. For 
the following we will detail the derivation for the 1 Sq state, but only quote the results for 
the 3 Si state. 

We first derive the fragmentation function Dq^p(z) with the leading terms in the HQET 
propagator and vertex only. The amplitude reduces to 



iM l = V s 1 ; ^- u(p') 1 + —i 7 5 (1 + i T . (12) 

3 m q [s — mq) \ n ■ k J 

Squaring and summing over spins and colors of the light quark, we get 

z(l-z) 



E| _ Mi|2 = 64^(0)1^, + 
9 mi v 



M 3 (l - (I - r)z) 2 {s - m 2 ) 2 



-l + z + 3rz ArM 



M{l-{l-r)z){s-mlf (s - m 2 Q )\ 
The corresponding amplitude-squared for producing an on-shell heavy quark is 



(13) 



£|M>| 2 = — Tr(rf(l+^)) . (14) 



z 



Substituting \J\A\ 2 and |A^o| 2 m to (P) and integrating over s, we get 

2a s (2m (? ) 2 | J R(0)| 2 rz 3 (l - zf 



8l7rm3 (l-(l-r)^) 6 
x (3(1 - (1 - r)z) 2 - 8rz(l - z) + Y2rz{\ - (1 - r)zj) . (15) 



Expressing this in terms of y using (|T]) and expanding to next-to-leading order in r, we have 



D Q ^ P {z) = Ar^-i)!^ 2 + Ay + 8) - N^^(3y 2 + Ay + 8) + 0(r) , (16) 



1 •:> i ■ 

ry v * " y v 



where A" = 2a 2 r |i?(0)| 2 /(8l7rm 3 J. Therefore, in terms of a(y) and b(y), the leading term in 
the HQET Lagrangian contributes 

a(y) = N^^(3y 2 + Ay + 8), (17) 
h(y) = N^-^(-(y-l)(3y 2 + Ay + 8)) . (18) 

The corresponding calculation for the 3 Si state gives a*(y) = 3a(y) and b\{y) = 3bi(y). 
These contributions to Dq^p(z) and Dq^v(z) differ by a spin factor of 3, as required by 
heavy-quark spin symmetry. 

Next we calculate the contributions from 1/mg corrections in the heavy-quark propagator 
and the heavy-quark vertex. Expanding out the I/ttiq correction to the propagator in ( [TTj) 
to first order, the correction to the amplitude is 



8y/^a s R(0) M 2 VM 1 / m q G 



m q (s - m l Q) 2 \ mq 2m c 
} ■ k 

7 

n ■ k 

Keeping the interference terms in \Mi + M 2 \ 2 , summing over spins and colors, and inserting 



x u( P ') { 1 + ^V) 7 5 (i + f) r . (19) 



into (§), we find a I/tuq correction to Dq^p(z). Expressing this in terms of y, we find that 
the contribution to b(y) is 

b 2 (y) = ((-2^ + C 2 y)(3y 2 + Ay + 8)) . (20) 

A similar calculation for the 3 Si state gives b* 2 {y) = 3b 2 (y). The 1/mg correction to the 
amplitude in (TTTJ) from the heavy-quark vertex is 



8J^a s \R(0)\M 2 VM 1 / n„k u \ , ,, .. 



Keeping the interference terms in \Aii + .M3I 2 , we obtain after some work the contribution 
to b(y) and to b*(y) due to the 1/mq vertex correction 

b 3 (y) = N^-fi- (-C 2 {y - \){3y 2 + Ay + 8) + 6d(y -l)(y + 2)- 12C 3 y) , (22) 

b* 3 (y) = 3N^-^- (-C 2 (y - l)(3y 2 + Ay + 8) + QC^y - \){y + 2) + 4tf 3 j,) . (23) 

In (p!l|), the (1 + ^)/2 factor adjacent to T projects onto the large component of the 

heavy-quark spinors produced by the source T. There is also a contribution of order I/tuq 

from the small component of the heavy-quark spinors of the fragmenting Q quark [[J . The 

corresponding amplitude is given by 

8^a s R(0) M 2 VM 1 / n^k u \ 
4 3 m q (s — nig) 2 n ■ k J 

x«(p')tV(i + f) (^Wt - fr")) ^y^ r • ( 24 ) 

The contributions to b(y) and to b*(y) from the interference term in \M.\ + M.^\ 2 are 

h(y) = 2N^^(3y 3 + 5y 2 + 2y-4), (25) 

= 6iV^l(y 3 - y 2 + 2y - 4) . (26) 

y 5 



The complete expression for b(y) is obtained by adding (|i8|) , (p0|) , (|22|) , and (p5|) . Thus 
the fragmentation function Dq^p(z) for the 1 5'o state, to next-to-leading order in I/thq, is 
given by @ with 

o(y) = N^^(3y 2 + Ay + 8), (27) 
6(y) = N^^((y-l)(3y 3 + 15y 2 + 8y-8)-12(C 3 -l)y 2 ) . (28) 



The complete expression for b*(y) is obtained by adding 3bi(y), 3b 2 (y), (|23"D , and (p6|) . The 
fragmentation function Dq^ v (z) for the 3 Si state, to next-to-leading order in I/tuq, is given 
by © with 

a*(y) = 3N^^(3y 2 + Ay + 8) , (29) 
b*(y) = 3N^-^(-(y-l)(y 3 + y 2 -8y + 8)+A(C 3 -l)y 2 ) . (30) 



The terms proportional to C 2 in (|28|) and (|30"D cancel between propagator and vertex cor- 
rections. We have set C\ = 1 in fl28|) and (|30D . If we further put C3 = 1, we recover the 
next-to-leading terms in the 1/r expansion of the PQCD fragmentation functions given in 
Ref. [§. 

The heavy-quark mass expansions (0) and (|3|) break down in the limit y — > 00, which 
corresponds to z — > 0, and also in the limit y — ► 1, which corresponds to z — > 1. As y — > 00, 
the leading terms, given by (E7J) and (59), scale like l/(ry 2 ), while the next-to- leading terms 
in (^) and fl30|) scale like l/y. Thus the 1/rriQ expansion breaks down when y is of order 
1/r or larger. In the limit y — > 1, the leading terms in (0) and ([3]) vanish like (y — l) 2 /r while 
the terms proportional to C3 — 1 in the next-to-leading terms go to as the first power of 
y — 1. Thus, unless C3 = 1, the expansion also breaks down for y — 1 of order r or smaller. 

In Fig. 0, we compare the PQCD fragmentation functions (solid curves) with the heavy- 
quark mass expansions (fj) and (|3|) at leading order (dotted curves) and next-to-leading 
order (dashed curves) in r. We use the value r = 0.10, which corresponds to D mesons. 
The normalization is fixed by arbitrarily setting iV = 1 in (PT|)-(P(1D. Note that we have set 
C3 = 1 in ( PB|) and (|SI1|). For any other value of C3, either (0) or ([5]) becomes negative for 
y very close to 1 indicating the breakdown of the 1/rriQ expansion when z is too close to 1. 
From the figure it is clear that the next-to-leading order curves are in very good agreement 
with the complete PQCD fragmentation functions for both D and D* mesons. Surprisingly, 
the leading order result for fragmentation into D* mesons also agrees very well with the 
complete PQCD result, while the leading order result for the D meson falls about 30% low 
near the peak. 

III. PQCD Model for Heavy-Quark Fragmentation 

It is tempting to use the heavy-quark limits of the PQCD fragmentation functions as 
phenomenological models for the fragmentation of a heavy quark Q into heavy-light mesons 
Qq, where Q = c or b and q = u,d, or s. To next-to- leading order in 1/rriQ, these fragmenta- 
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tion functions are given by @ and (|3|), with a(y), b(y), a*(y), and b*(y) given in ( p7[) - (|30|). 
In addition to iV and r, we must treat C3 as a phenomenological parameter, since, according 
to (^|), it depends on the low-energy scale /1 where perturbation theory breaks down. These 
3 parameters all have well-defined scaling behavior with the heavy-quark mass; namely, N 
is independent of mg, r scales like 1/mq, and C3 scales like a s (mQ) 9 ^ 33 ~ 2n f\ Thus, if the 
parameters are determined phenomenologically from data on charm fragmentation into D 
and D* mesons, then the corresponding parameters for the B and B* mesons can be deter- 
mined by scaling. The problem with this model is that unless C3 = 1, either Dq^p(z) or 
Dq^v(z) becomes negative for z near 1. This unphysical behavior only arises in a region of z 
where the 1 / tjiq expansion is breaking down, but it makes these fragmentation functions less 
attractive as a phenomenological model. If we choose C3 = 1 to avoid these difficulties, we 
might as well avoid the 1/rriQ expansion altogether and use the complete PQCD fragmenta- 
tion functions as our model. We therefore propose as a model of heavy quark fragmentation 
the PQCD fragmentation functions calculated in Ref. B: 



Tzi\ — zS 2 r 

D Q ^p(z) = N - ) J -— 6-18(l-2r)z+(21-74r + 68r 2 )z : 



v l - (1 - r)zf 

-2(l-r)(6-19r + 18r 2 )2 3 + 3(l-r) 2 (l-2r + 2r 2 )z 4 ] , (31) 
Dq->v{z) = SN _ l (i _ J [2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )z 2 

-2(l-r)(4-r + 2r 2 )^ 3 + (l-r) 2 (3-2r + 2r 2 )z 4 ] . (32) 

The only parameters are the normalization N, which is independent of ttlq, and r, which 
scales like l/rriQ. The parameter r, which in the PQCD calculation has the value m q /(mQ + 
m q ), can be interpreted as the ratio of the constituent mass of the light quark to the mass 
of the meson. Integrating over z, we obtain the total fragmentation probabilities: 
"i /« + 13r + 228r 2 -212r 3 + 53r 4 



/ dzD Q ^ P {z) = 3N 
Jo 



15(1 -r) 5 
r(l + 8r + r 2 — 6r 3 + 2r 4 ) log(r) 



(1 -r) 6 

"i /24 + 109r - 126r 2 + 174r 3 + 89r 4 



(33) 



/ dzD Q ^ v {z) = 3N 
Jo 



15(1 -r)' 
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r(7 - 4r + 3r 2 + 10r 3 + 2r 4 ) log(r) 



+ ^'~— J- >»™> (34) 



The PQCD fragmentation functions ([31]) and fl32|) give the distributions in the longitu- 
dinal momentum fraction z for the mesons P and V in a heavy-quark jet. This model can 
be easily extended to give the distribution in their transverse momentum kp relative to the 



jet momentum |^T|. In Ref. |J, the fragmentation functions were obtained as integrals over 
the invariant mass s of the fragmenting heavy quark: 

ds 



D Q ^ P/V (z) = — d Q ^ P/v (z, s) , (35) 

where the lower limit of the integration is: 



, , M 2 r 2 M 2 , . 

Smin(^) = h : _ • (36) 

Z X z 



The functions dQ^p/y(z, s) in the integrand are given by 

(1 -z)(l + rz) 2 



d Q ^ P (z,s) = 6NM z rs 



(1 - (1 - r)z) 2 (s - (1 - r) 2 M 2 ) 2 
[2(1 - 2r) - (3 - 4r + Ar 2 )z + (1 - r)(l - 2r)z 2 ]M 2 4r(l - r)M 4 



d Q ^ v (z, s) = 6NM 2 rs 



(1 - (1 - r)z)(s - (1 - r) 2 M 2 f (s - (1 - r) 2 M 2 ) 4 _ 

(37) 

(1 - z){l + 2rz + {2 + r 2 )z 2 ) 



(l-(l-r)z) 2 (s-(l-r) 2 M 2 ) 2 
[2(1 + 2r) - (1 + 12r - Ar 2 )z - (1 - r)(l + 2r)^ 2 ]M 2 12r(l - r)M 4 



(1 - (1 - r)z)(s - (1 - r) 2 M 2 ) 3 (s - (1 - r) 2 M 2 ) 4 



The invariant mass s is related to kp and z by 



2 1 — z 

where M = mQ+m q in the nonrelativistic limit. If, instead of integrating over s, we integrate 
over z with k\ held fixed, we obtained the kp distribution for the fragmentation process. 
Introducing the dimensionless variable t = kp/M, we can define the /^-dependent functions 
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d Q ^p/v{z,t) and D Q ^ P/V (t) by 

dt D Q ^ P / V (i) = dt dz d Q ^p/ v (z,t) 



1 roo dg 

dz — d Q ^ P/v (z, s) . 

J s min( z ) & 



This implies 



D Q ^P/v(t) = 2M 2 t f 1 dz 

J 



o z(l — z)s(z, t) 



d Q -> P / v {z,s(z,t)) 



with 



s(z,t) = M 2 



1+t 2 r 2 + t 2 



1-z 



Carrying out the integrals over z, we find 



Dr 



,(t) 



24rt 



4r - (2 + r + 2r 2 )t 2 log(r 



(40) 
(41) 



(42) 



Nr 1 
2(1 -r) 6 ¥ 

(1 - r)t [30r 3 - r(61 - 20r + 28r 2 )t 2 - (3 - 48r + 48r 2 - 12r 3 )t 4 



+ Ut [4r 3 - r(2 + r + 2r 2 )t 2 + (1 - r) 2 t 6 ] log 
+ 3 



V 2 + t r 



, 1 + t2 / 

10r 4 - 3r 2 (ll + 2r + 2r 2 )t 2 + (3 + 4r + 19r 2 - 6r 3 )t 4 



-(3 + 12r - 20r 2 + 8r 3 )t 6 



Arctan 



(1 -r)t x 
r + t 2 



(43) 



3Nr 



-8rt 



12r 2 - (6 + 7r + 2r 2 )t 2 log(r 



2(1 -r) 6 t 6 

(1 - r)t [30r 3 - r(61 + 28r - 20r 2 )t 2 + (5 - 8r + 8r 2 + 4r 3 )t 4 



+ 4t 
+ 



12r d 



r(6 + 7r + 2r 2 )r + (1 - r) 2 r log 



V 2 + t 2 * 



1 + t 2 J 

30r 4 - 3r 2 (33 + 22r - 10r 2 )t 2 + (9 + 20r + r 2 + 22r 3 + 8r 4 )t 4 



+ (9- 12r + 4r 2 + 8r 3 )£ 6 



Arctan 



(l-r)^ 
r + t 2 



(44) 



In general, fragmentation functions D(z,fi 2 ) depend not only on z but also on a factor- 
ization scale /i. In a high energy process that produces a jet with transverse momentum pt, 
the scale /i should be chosen to be on the order of pt- The functions fl3~l|) and (|32|) should 
be regarded as models for heavy-quark fragmentation functions at a scale /i of order tuq. 
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For values of \x much larger than mq, the fragmentation functions fl3T|) and (|32|) should be 

evolved from the scale tuq to the scale \i using the Altarelli-Parisi equation: 

d dvt z 

a 2 — D Q ^ H (z, /i 2 ) = J z J Pq^q(~, //) D Q ^ H (y, /i 2 ) , (45) 

where Pq^q(x) is the appropriate splitting function: 

„ . . 2a s (u) (l + x 2 \ , . 

W^) = ^(— J + . (46) 

One aspect of the initial conditions ( p?l"D and (p?2| ) and the evolution equation (KSj) that 
may cause problems in practical applications is that they do not respect the phase space 
constraint: 

D Q ^ H (z^ 2 ) = for z<M7/i 2 , (47) 

This can be remedied || by using (47) as the initial condition on the fragmentation function 

equation and replacing ([45]) by the inhomogeneous evolution equation 

d dy z 

H 2 -^Dq^h(z,h 2 ) = / — Pq^ q (-,ii) DQ^ H (y,yjj, 2 )+d Q ^ H (z,n 2 )e(jj 2 -s min (z)) , (48) 
Ofj, J z y y 

where dQ^n( z , s) is defined by the integrand in ( p5|) and s min (z) is given in (|3 



IV. Comparison with other Fragmentation Models 

The model for heavy-quark fragmentation which has been used most extensively in phe- 
nomenological applications is the Peterson fragmentation function [|lC|j : 

where Nh and e# are adjustable parameters that may depend on the hadron H. This 
fragmentation function has the correct behavior in the heavy-quark limit if Nh scales like 
1/mq and e# scales like 1/rriQ. Identifying e# with r 2 and expressing (|4*9| ) in terms of the 
Jaffe-Randall scaling variable y defined in ([!]), we find that it reduces in the limit r — > to 

D Q->»M ~> — [( y -l)2 + l]2 " ( 50 ) 
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The Peterson fragmentation function is just the square of a light-cone energy denominator 
multiplied by a phase space factor. It contains no spin information; the normalization 
parameter Nu is to be determined independently for the pseudoscalar and vector mesons of 
a heavy-quark spin multiplet. 

An alternative fragmentation model which does contain spin information has been pro- 
posed by Suzuki jnj. Suzuki's fragmentation functions are derived from the same Feynman 
diagram in Fig. [I] as the PQCD fragmentation functions, but with two essential differences. 
First, the diagram was calculated in Feynman gauge. If a general covariant gauge had been 
used, Suzuki's fragmentation functions would have depended on the gauge parameter. The 
PQCD fragmentation functions that we calculated are gauge-invariant. We calculated the 
diagram in the axial gauge only for simplicity. If we had used a covariant gauge, we would 
have had to also include diagrams in which both the virtual heavy quark and the virtual 
gluon are emitted by the source V in Fig. |I[ Alternatively, we could have calculated the 
PQCD fragmentation functions for the fragmentation of a heavy quark into S-wave heavy 
quarkonium directly from the general gauge-invariant definition ||12|| . Such a calculation has 
been carried out for the equal mass case of charmonium by Ma [IB]. A second essential 



difference between the PQCD model and Suzuki's is that we integrated over the invariant 
mass s of the fragmenting quark (see Eq. (^)). The invariant mass is related to the transverse 
momentum kx of the meson relative to the fragmenting quark by (^). Rather than integrat- 
ing over k?, Suzuki chose to evaluate the integrand at a typical value (kj)- Suzuki's model 
therefore has 3 parameters: the overall normalization N, the mass ratio r, and (k^/rriQ. 
When expressed in terms of the scaling variable y defined in (HI), Suzuki's fragmentation 
function Dq^p(z) reduces in the limit r — > to 

N ^- 2 ) 2 + k 2 

Aj-p(z) - - (y - 1) [y2 + ft2]4 » ( 51 ) 

where k 2 = (k^) / (r^) . By heavy-quark spin symmetry Dq^y{z) differs, in this limit, 
only by a factor of 3. 

The Peterson, Suzuki, and PQCD fragmentation functions all vanish like (1— z) 2 as z — ► 1. 
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An alternative fragmentation function which vanishes like the first power of (1 — z) has been 



proposed by Collins and Spillers [14]. This was motivated by incorrect dimensional counting 



rules. The correct dimensional counting rules for QCD [15] do in fact give a limiting behavior 
of (1 — z) 2 for the fragmentation function. The Collins- Spillers fragmentation function can 
be derived in a similar way to ours, except that in the Feynman diagram in Fig. ffl, the 
interaction mediated by the virtual gluon is replaced by a point-like scalar Yukawa coupling 
between the meson, the heavy quark, and the light quark. Consequently, the denominator of 
the matrix element contains only one power of (s — Trin), in contrast to the 2 powers in (fl2|). 
It is the omission of the gluon propagator that changes the behavior as z — > 1 from (1 — z) 2 
to (1 — z). Also, instead of integrating over the invariant mass of the fragmenting quark as in 
(0), Collins and Spillers, like Suzuki, evaluated the integral at a typical value {kj.). Taking 
the scaling limit r — > 0, the fragmentation function of Collins and Spillers reduces to 

Dq^p(z) - — (V " 1) {y [y2 I k2]2 , (52) 

where k 2 = (fcf.) / (r 2 mg) . 



The various fragmentation models in the literature have been summarized in Ref. |16 
and compared with experimental data on D and D* production. The string models and 
parton cluster models are very different in spirit from those discussed above. One can derive 
analytic expressions for the heavy-quark fragmentation functions from the string models [17|. 
They contain a tunneling factor exp(— Bm 2 H /z), which suppresses the small- z region. In the 
scaling limit, the Lund symmetric fragmentation function behaves like 

D Q ^ H (z) -> Nr^e~ B{m ^ +{k ^ (y - if . (53) 

Unless N scales like e Bm «mg +1 , this is inconsistent with heavy-quark symmetry, which re- 
quires the leading term to scale like mq as mq — > oo. 

The PQCD model for heavy-quark fragmentation has a number of advantages over those 
described above. First, it is rigorously correct in the limit m q ^> Aqcd- Higher order 
perturbative corrections can be systematically calculated. Relativistic corrections can also 
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be calculated in terms of additional nonperturbative matrix elements [18]. Second, our model 
is consistent with heavy-quark symmetry in the limit mg — > oo. The logarithms of mq that 
are predicted by HQET would be reproduced by the higher order perturbative corrections. 
The PQCD model is also more predictive than those in Refs. 11, [14[]. It describes spin- 
dependent effects, like Suzuki's model, but without introducing any additional parameters. 
The PQCD model not only predicts the ^-dependence of the fragmentation functions but 
also their dependence on k?, the transverse momentum of the meson relative to the jet. The 
fragmentation functions ([31]) and ([32]) apply only to S-wave mesons, but the fragmentation 
functions for higher orbital-angular-momentum states can also be calculated. The PQCD 
fragmentation functions for the P-wave mesons have been calculated to leading order in a s 



in Refs. 19 



V. The Vector-to-Pseudoscalar Ratio 

In any production process for heavy-light mesons, one of the most fundamental experi- 
mental observables is the ratio 

Py = VTp ' (54 > 

which measures the relative number of vector mesons V and pseudoscalar mesons P that are 
produced. If the mesons are produced within a heavy-quark jet, then V and P in fl54[) can be 
identified as the fragmentation probabilities for the heavy quark to fragment into vector and 
pseudoscalar mesons, respectively. The ratio Py can depend on kinematic variables, such as 
the longitudinal momentum fraction z of the meson or its transverse momentum fc^ relative 
to the jet. In the PQCD model for fragmentation, the normalization factor iV cancels out 
in the ratio (|54l) , so that Py is determined by the parameter r only. 

The simplest measure of the ratio Py comes from the total numbers of vector and pseu- 
doscalar mesons in the jet integrated over z and kx- Setting P and V in fl54[) to the frag- 
mentation probabilities in fl3~3|) and (|34]) , we find that the ratio Py in the PQCD model of 
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fragmentation is 

;i - r)(24 + 109r - 126r 2 + 174r 3 + 89r 4 ) + 15r(7 - 4r + 3r 2 + 10r 3 + 2r 4 ) log(r) 



Pi 



2(1 - r)(16 + 61r + 51r 2 - 19r 3 + 71r 4 ) + 60r(2 + r + r 2 + r 3 + r 4 ) log(r) 

(55) 

This ratio is plotted as a function of r in Fig. |3|. From the graph it is clear that Py is not 
strongly dependent on r. At r = 0, Py = 3/4 as required by heavy-quark spin symmetry. 
As r increases Py decreases slowly to Py = 0.51 at r = 0.5. Thus at nonzero values of 
r, the vector state is less populated than would be given by naive spin counting. We can 
determine the value of r for the D and D* system using experimental measurements of Py. 
A complete compilation of experimental data for Py from LEP, CLEO, ARGUS, PETRA, 
and TRISTAN can be found in Ref. ]20|]. The key point in obtaining consistency between 



these measurements is using the updated branching ratio B(D + * — ► D tt + ) « 0.68 instead of 
the old value 0.55. The experimental value Py = 0.65 ± 0.06 determines the parameter 
for the D — D* system to be ro = 0.10 _ao7- If we interpret r as the ratio of the constituent 
mass of the light quark to the mass of the meson, then the value rp = 0.10 corresponds to 
a constituent mass of 200 MeV. Given a value of rp, we can determine the corresponding 
value for the B — B* system by using the simple scaling behavior vb = (mo/mB^D- This 

+0.04 

gives r B = 0.03 _ .02- 

Having determined the parameter r from data on D — D* production, we can now predict 
how the vector-to-pseudoscalar ratio should vary as a function of the longitudinal momentum 
fraction z. The ^-dependent ratio Py(z) is defined by ([54]), where P and V are given by the 
fragmentation functions (|3l|) and (p2|) : 

with 

n(z) = 2-2(3-2r)z + 3(3-2r + 4r 2 )^ 2 

-2(1 - r)(4 - r + 2r 2 )^ 3 + (1 - r) 2 (3 - 2r + 2r 2 )z 4 , (57) 
d(z) = 3-3(3-4r)z+ (12 - 23r + 26r 2 )z 2 



-(1 - r)(9 - llr + 12r 2 )z s + 3(1 - r) 2 (l - r + r 2 )z* . (58) 

This ratio is plotted as a function of z in Fig. [| for the values r = 0.10 (solid curve), r = 0.03 
(dotted curve), and r = 0.22 (dashed curve). At z = 0, -FV^) = 1/2, regardless of the value 
of r. It decreases slightly for small z, and then increases monotonically to a maximum value, 
at z = 1, of 0.74 for r = 0.03, 0.73 for r = 0.10, and 0.70 for r = 0.22. Note that, in spite of 
the large uncertainty in our determination of r, the uncertainty in Py(z) is less than about 
11%. Thus the PQCD model gives a rather unambiguous prediction that Py should vary 
from around 1 /2 at small values of z to almost 3/4 near z — 1 . 

The ^-dependent ratio Py{kr) is defined by (151), with P and V given by (f4"3"D and 



= 3 wi + n 2 log(r) + n 3 log (gg) + n 4 Arctan (j^g) 
y T 4 di + ^2 log(r) + d 3 log (§f ) + rf 4 Arctan ' 



where 



and 



ni = -(l-r)t[30r 3 -r(61 + 28r-20r 2 )t 2 + (5-8r + 8r 2 + 4r 3 )t 4 ] , (60) 

n 2 = -8rt[12r 2 - (6 + 7r + 2r 2 )t 2 ] , (61) 

n 3 = 4t[12r 3 - r(6 + 7r + 2r 2 )t 2 + (1 - r)H 6 } , (62) 
n 4 = [30r 4 -3r 2 (33 + 22r-10r 2 )t 2 + (9 + 20r + r 2 + 22r 3 + 8r 4 )t 4 

+ (9-12r + 4r 2 + 8r 3 )t 6 ] , (63) 

rfx = -(1 - r)t[30r 3 - r(61 + 16r - 8r 2 )t 2 + 3(1 + 2r - 2r 2 + 2r 3 )t 4 ] , (64) 
d 2 = -12rt[8r 2 - 2(2 + 2r + r 2 )t 2 ] , (65) 
d 3 = 6t[8r 3 -2r(2 + 2r + r 2 )t 2 + (l-r) 2 t 6 ] , (66) 
d A = 30r 4 -9r 2 (ll + 6r-2r 2 )t 2 + 3(l + r) 2 (3 + 2r 2 )t 4 + 3(3-4r 2 + 4r 3 )t 6 . 

(67) 

This ratio is plotted as a function of t = krp/M in Fig. |5| for the three values r = 0.10, 0.03, 
and 0.22. At t — 0, Py(t) = 3/4, regardless of the value of r. As t increases, Py(t) quickly 
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decreases to its asymptotic value at t = oo. At t = 1, Py(£) is within 0.1% of its asymptotic 
value of 0.65 for r = 0.03, 0.62 for r = 0.10, and 0.60 for r = 0.22. Again we find that, in 
spite of the large uncertainty in r, we obtain a rather precise prediction for Py as a function 
of krp- 

The PQCD fragmentation functions for vector mesons have been applied previously PI1 
as a phenomenological model to describe the fragmentation processes c — > D* and b — > B* . 
The fragmentation functions were separated into the transverse and longitudinal polarization 
components. The spin alignment, which measures the ratio of transverse to longitudinal po- 
larizations, was calculated as a function of z and as a function of kx- In the case of production 
of D* by charm fragmentation, the spin alignment predicted by the PQCD fragmentation 
model was shown to be consistent with CLEO measurements [^TJ. In addition, the predicted 
value of the average longitudinal momentum fraction (z) for c — ■> D* and for b — > B* was 



shown to be in excellent agreement with data from LEP, CLEO, and ARGUS [21 1. The 
values of r used for D* and B* mesons in these comparisons were r = 0.17 and r = 0.058, 
respectively, which lie within the range determined above from measurements of Py. 

The PQCD fragmentation functions have also been applied in Ref. [P^l to predict the 
fragmentation spectra for the B s and B* mesons based on the production rates of the B s 
mesons measured at LEP. Instead of treating the normalization N as a phenomenological 
parameter as advocated in this paper, the authors calculated N using the PQCD expression, 
which involves a s at the scale of the strange quark mass. 



VI. Summary 

We have studied the heavy-quark mass limit of the PQCD fragmentation functions for 
producing S-wave mesons. The leading and next-to- leading terms in I/tuq were calculated 
directly from HQET. The PQCD fragmentation functions were proposed as a phenomeno- 
logical model for fragmentation into heavy-light mesons. With only 2 parameters, this model 
describes fragmentation into the 1 5'o pseudoscalar meson state and the transverse and longi- 
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tudinal polarization states of the 3 5i vector meson. It describes not only the z-dependence 
of the fragmentation probabilities, but also their dependence on the transverse momentum 
kr of the meson relative to the jet within which it is produced. This model can easily 
be extended to describe heavy quark fragmentation into P-wave states using the PQCD 
fragmentation functions calculated in [Tj|. The PQCD fragmentation functions were com- 
pared with other models for heavy-quark fragmentation in the literature. As an application, 
the PQCD fragmentation functions were used to predict the ratio of vector-to-pseudoscalar 
states as a function of z and as a function of kx- The ratio Py is predicted to vary from 
around 1/2 at small values of z to almost 3/4 near z = 1. 
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Figure Captions 

1. Feynman diagram used to calculate the PQCD fragmentation functions in axial gauge. 

2. Comparison of the D c _ D (z) (lower set of curves) and D c ^ D *{z) (upper set of curves) 
fragmentation functions. The normalization is arbitrary. Shown are the full PQCD 
results (solid curves), the leading terms (dotted curves) in the heavy-quark mass ex- 
pansion, and the leading plus next-to-leading terms (dashed curves) in the heavy-quark 
mass expansion. 

3. The ratio P v as a function of r. 

4. Predictions for the ratio Pv(z) as a function of z for r = 0.10 (solid curve), r = 0.03 
(dotted curve), and r = 0.22 (dashed curve). 

5. Predictions for the ratio Pv(kr) as a function of kr for r = 0.10 (solid curve), r = 0.03 
(dotted curve), and r = 0.22 (dashed curve). 



24 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409316vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409316vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409316vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409316vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409316vl 



